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Background 


* Cold pool activity can significantly impact horizontal wind velocities at 10-30 
km grid scale, leading to significant Inhomogenities in GCMs. (Tompkins 
2001; Ross et al. 2004; Zuidema et al. 2012). 


¢ Elsaesser et al. (2013) utilized the QuikSCAT data over tropical oceans to 
assess the relationship between the surface winds and cold pools using cold 
lelele) mM ebeleinlomoelantaia(O) GaP 


¢ CPKE used the anomalies in the horizontal wind and since variance in surface 
winds over tropical oceans is likely related to cold pool outflow or mesoscale 
downdrafts from nearby regions of convection. 


¢ Elsaesser et al. (2013) did not 
consider tain contamination and 
assumed that CPKE was capturing 
7Xoiuee-limvishale mn hal-inleselsme-bxelenalemmnele 
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* Kilpatrick et al. (2015) used _ Sutace Wind 
ASCAT data after limiting rain ‘Stee 
contamination to observe MCS : = 
downdrafts using divergence. = a aa 


Cold Pool Kinetic Energy 


° Kilpatrick et al. (2015) found that . A ie 
ASCAT cannot detect the shusetrieisioc cemsonnsin 
convective-scale gust fronts, but & i 
can detect mesoscale downdrafts 
corresponding to cold pools. 
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Courtesy: Elsaesser et al. (2013) 


Motivation and Goals of this study 


¢ We hypothesize that an approach based on gradient features (GEFs) can 
be used to identify the cold pools over tropical oceans. 


* Cold pools form a boundary of gust front, thus creating an area of steep 
oradients in horizontal winds. 


¢ Divergence and vorticity are noisier to identify the areas of steep 
oradients in the horizontal wind, thus gradient of horizontal wind is 
calculated using: 


Ou n OV 
Vi = ihe Ose 
~~ 10u OV 


dy dy 


¢ After calculating gradient in the horizontal wind, Sobel technique and 
convex hull algorithm was implemented to identify the perimeters of the 


GEFs. 


¢ Properties corresponding to individual Gs like eccentricity, area, 
weighted centroid etc. were calculated. 


¢ Identification of the boundary of the density current was performed 
using divergence (V.V) and vorticity VXV). 


¢ WRE simulation was ran and Gs were obtained using model winds to 
test our hypotheses. 


Data used in this study 


Observations 


¢ ASCAT Metop-A dataset from Remote sensing systems (REMSS) its used at 12.5km resolution 
for the period of 2007-2015. 


¢ This dataset uses TRMM Microwave Imager (IMI), GPM Microwave imager (GMI), Advanced 
Microwave Scanning Radiometer-2 (AMSR-2) and WindSat to determine if rain is present within 
GW aloleracw- lem elem (olertaleseMos aval a0 Mol oarisloser 


* CPC Morphing Technique (CMORPH) precipitation dataset for the above period at 8km and 30- 


min resolution is used. 


¢ Modern-Era Retrospective analysis for Research and Applications (MEERA-2) products were 
used to identify background fields. 


e S-Polka Radar data from DYNAMO dataset archive was used to visualize convection over Indian 
Ocean corresponding to ASCAT swaths and WRF simulations. 


Modeling Framework 


¢ WRE with ARW core was run on a nested grid of 27-9-3 km using initial 
conditions from GFS for the period of Active MJO-1 (19% and 20% 
October 2011) with following schemes: 


Longwave Radiation RRTMG 

Shortwave Radiation RRTMG 

Microphysics Thompson 

Boundary Layer YSU 

Cumulus Tiedtke (27km and 9km) 


* Convection was allowed explicitly in 3km domain. 


NY Card skeye Co) Coyeay 


¢ Post the identification of Gs globally, a comparison was drawn between 
ASCAT and MERRA-2 to see whether the reanalysis dataset is able to 
resolve the background dynamics. 


¢ Surface sensible and latent heat fluxes were calculated using MERRA-2 
|oy-Tol <cesnoyebaels s(l(em oyaolelelolncw-bele mate Ova MR vsnelely 


° Precipitation from CMORPH was collocated with the gradient features 
within 8 degrees of the centroid. 


¢ WRE’s 3 km dataset was regridded to 12.5 km and GEFs were calculated 
for the domain. 


¢ To study the thermodynamic structure of GEFs, virtual temperature (T,) 
was calculated from model parameters and GI's were overlaid on T, to see 
whether they collocate or not. 


Results 
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Fig.1 Density of GFs in a 0.5° grid box for 2007-2015 
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GFs in the Indian Subcontinent (19-20 October 2011) 


Fig. 3 GFs over the Indian Subcontinent for 19 — 20 October 2011 (DYNAMO MJO-1) 
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Fig.4 Divergence and Vorticity for ASCAT (left) and MERRA-2 (right) 
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Fig.5 Surface Latent heat and sensible heat flux (Wm7~) for ASCAT (left) and MERRA-2 (right) and CMORPH 
precipitation (mm/hr) 
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Fig.6 Reflectivity (dBZ) and Relative velocity (ms*) for 20 October 2011 16:32:21 UTC as 
eater ibinexemineseamaXerchu mace) melutunetoa OD ANv.WULOM@AIIED 
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Fig.7 Divergence, vorticity and fluxes from ASCAT and MERRA 


Fig.8 GFs over WRF domain for 19™ -20" October 2011 
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Frequency of Gradient features in Indian Ocean from WRF 
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Fig. 9 Frequency of GFs over WRF domain from ASCAT (left) and WRF (right) for 
19 -20" October 2011 
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¢ For 2007-2015, global density of Gs shows that these features match well 
with the precipitation pattern over the tropical oceans, which gives us a 
confidence in our hypothesis. 


¢ Divergence and vorticity corresponding to the boundary of observed GFs 
gives a good picture of the outflow boundary of a mesoscale system. 


¢ Surface sensible and latent heat fluxes calculated using ASCAT winds produce 
a better result as compared to MERRA-2 inherent fluxes. 


¢ Precipitation around the perimeter of the features shows that these features 
are forming along the wake of the density current. 


* Reeridded WRF is creating smaller cold pools overall, as compared to the 
climatological distribution of GEFs, which might have been due to steeper 
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¢ WRF is not able to resolve the observed GFs very well, which might also 
be another reason that its producing smaller GFs. 


¢ Low Virtual Temperature (Tv) and GF are collocating well in WRF, thus 
suggesting that these features are rather cold pools. 


Future Work 


¢ We are working on a tool to match the radar and ASCAT swath for 
DYNAMO, so as to create a better picture of convection in the region. 


° We will be running WRF with ECMWF initial conditions and will be 
comparing the results with the current simulation. 


¢ We will also be including the winds from ASCAT and RapidSCAT in 
WRE and will study the difference in the simulations. 


¢ We will be including more thermodynamic and dynamic parameters to 
study the cold pools over the global oceans. 


¢ We will also try to run our algorithm using a mesoscale model to see 
whether it can resolve the cold pools comparable with ASCAT. 
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